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ABSTRACT 

This computer simulation investigates the ways in which an embedding 
lattice influences a collision event. Low energy events were studied with 
a Born-Mayer potential function; high energy with a Thomas-Fermi-Firsov 
potential. Conclusions were: (1) the lattice will increase ranges for 
particle energies above 400 ev; (2) the apparent mass concept is not a 
valid description of events in a lattice; (3) lattice effects will signi- 
ficantly modify the low energy portion of the target atom energy distri- 
bution function; (4) there is no evidence that a copper atom will "rebound" 
from a copper lattice, 

The writer wishes to express his appreciation for the assistance and 
encouragement given him by Associate Professor Don E. Harrison, Jr. of the 


U. S. Naval Postgraduate School in this investigation. 
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1. History. 

Radiation damage is the overall term applied to erosion, disruption, or 
rearrangement of a crystal lattice by bombardment of ions or atoms. Sput- 
tering, a phenomena sometimes observed in radiation damage, applies to the 
process of lattice atoms being ejected or "knocked out" of the crystal. The 
sputtering process was first observed in 1852 by Grove? in the cathodes of 
electrical discharge tubes, 

Since Grove's observation, experimental and theoretical physicists have 
been attempting to develop a sputtering theory that will predict measurable 
experimental quantities as a function of mass ratio, position in the periodic 
table, energy of bombarding particle and angle of incidence. To date, no 
such complete theory is available. 

A tremendous volume of material has been written on the subject of sput- 
tering since its discovery. [In 1955, Ervine listed 306 papers concerned 
with the sputtering process. These early papers disagree violently because 
important experimental conditions were not obtained or simply ignored; 
consequently many of the conclusions reached are contradictory. 

During this period two important theories were proposed. Kingdon and 
eee in 1923, developed a momentum transfer theory which agreed quite 
well with experimental data they had obtained for thoriated tungsten fila- 
ments. The "evaporation theory" was proposed in 1928 by von Hippel and 


4 
Blechschmidt and improved by Wore in 1944, It assumes that the incident 


. R. Grove, Phil. Trans. Roy. Soc. London, 142, 87 (1842). 
. M. Irvine, Dissertation, Lehigh University (1955). 


- von Hippel and E. Blechschmidt, Ann. Physik 86, 1006 (1928). 
- H. Townes, Phys. Rev. 65, 319 (1944). 
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particles heat a small volume of the crystal to a very high temperature 

and subsequently surface atoms evaporate. More recently (1953), Harrison’, 
following a suggestion made by Keywell’, presented a statistical theory of 
sputtering which could be placed somewhere between the two models mentioned 
above. We will not present details of these theories. We must mention how- 
ever, that they do not provide an adequate description of the sputtering 
process. 

The trend in sputtering theories has recently been toward the considera- 
tion of the individual collision processes within the lattice. Henschke® 
has used this approach in his theory of sputtering and Harrison and Mester: 
have also applied it, with modifications, to the theoretical study of sput- 
tering threshold energies. 

The theoretical work performed by Henschke uses the momentum transfer 
concept initiated by Kingdon and Lane to account for all the experi- 
mental phenomena observed at low incident particle energy. His basic treat= 
ment of collisions between atoms in a lattice is based on the assumption 

..ethat the collisions described can be treated with the general 

principles of classical mechanics, using impulsive forces, in a 

manner similar to the well-known collisions with restitution. 

However, he did not in all cases use the masses of the individual particles 
involved in the equations pertaining to the two body collisions. Instead 


he postulated an "effective" mass. If the collision is between a surface 


D. E. Harrison, Jr., Phys. Rev. 102, 1473 (1956). 

F. Keywell, Phys. Rev. 87, 160 (1952). 

E. B. Henschke, Phys. Rev. 106, 737 (1957). 

D. E. Harrison and G, D. Magnuson, Phys. Rev. 122, 1421 (1961) 
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atom and a moving ion or atom and is such that the surface atom is struck 

on its "inside" hemisphere, according to Henschke, an effective mass is not 
required and the two particles involved can be considered to have their actual 
masses, Should the collision of an atom or ion be directed inward from the 


target lattice surface, Henschke states, 


The bulk of the target is behind the struck atom and produces a 
very large ‘effective' mass compared to the mass of the ion or 
to the mass of the target atom,!1 


The reason Henschke did not assume pure elastic collisons was stated as 


Energy losses are due to the fact that the target atom is coupled 
rather strongly with the atoms of the lattice. Before the moment 
of highest compression is reached, the ion and the struck atom 
exchange energy with the neighboring atoms of the lattice. Debye 
waves are thus excited and dissipated irreversibly into the lat- 
tice, 


He goes on to state, 


..ethe final step in each sputtering process at any angle of inci- 
dence of the ion can be generally described as a collision of the 
ion with an upper surface atom, in which this atom is hit on its 
inside hemisphere so as to obtain an impulse with a component in 
the direction of the outward normal to the surface. If the energy 
transferred in this direction to this target atom by the impact 

of the ion is equal to or greater than the heat of vaporization, 
with which the atom is assumed to be bound to the crystal lattice 
plane, then this atom is ejected in the collision. 


Early study in radiation damage was precipited by the advent of the 
nuclear reactor. Recently however, with the imminent possibility of thermo- 


nuclear power production and ion propulsion engines for space vehicles, other 


practical applications of the processes involved in radiation damage have 


ie B. Henschke, Phys. Rev. 121, 1290 (1961). 
a oee reference 8, p. 738. 


oes reference 8, p. 737. 





become extremely important. 

Robinson and his coeworkenal” have made theoretical studies of the 
ranges in solids of atoms having energies from 1 to 100 kev using digital 
computer techniques on the basic assumption that the moving atom loses its 
energy through repeated binary elastic collisions with atoms of the solid. 
The masses used in these calculations are the true masses of the interact-~- 
ing particles, but they note that their assumption is certainly not valid 
below 100 ev. Gibson et. ee have also taken advantage of the speed avail- 
able with modern digital computers in the study of radiation damage. They 
do not assume binary collisions but instead use iteration techniques employ- 
ing Newton's equation of motion to solve the complex many body problem, The 
initial success of these programs have encouraged their originators to ex- 
plore further the possibilities of these techniques, Although the computer 
programs are designed to study radiation damage, the basic principles involv- 
ed are also important in sputtering. 

The work just discussed is very sensitive to the mathematical form of 
the interatomic potential. The Born-Mayer, Bohr (screened Coulomb), and 
the Thomas-Fermi-Firsov are some of the potentials used in these calcula- 
tions. The ability of any of the potential functions to describe a physical 
situation is highly dependent upon the energy range under consideration and 
the specific atoms or ions involved. No single potential has yet been de- 


vised that satisfactorily represents the interaction under all circumstances, 


peor S. Oen, D. K. Holmes, and M. T. Robinson, Jour Appl. Phys. 34, 
302 (1963). 


aa B. Gibson, A. N. Goland, M. Milgram, and G. H. Vineyard, Phys. 
Rev. 120, 1229 (1960) 
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An obvious unresolved question is inherent in these theories. The 
many body approach and the two body assumption, with a minimum energy limita- 
tion, appear to be reasonable approaches to the problem. Both methods des- 
cribe possible events in a lattice, but are they compatable? The general 
use of the two body approach is desirable because of its simplicity. The 
problem to be studied is: when can the two bedy collision be assumed, if 


at all, and should an "apparent" mass, as proposed by Handehke’, be employed? 





2. Objective. 

We intend to investigate theoretically the energy transfer process 
between an incoming copper atom and a copper lattice. The many body ap- 
proach, using computer techniques similar to Gibson et. as will be 
utilized and then compared with the simple two body solution. We hope 
the results will provide some quantitative answers to the questions posed 


in the preceding section, 





3. Procedure, 

The representative crystal lattice consists of a cube of 63 atoms 
arranged in a face centered cubic structure. This arrangement gives a 
total of 63 atoms. These atoms are free to move, if disturbed. [In an 
effort to simulate a larger lattice, the movable 63 atoms are surrounded 
on all sides by stationary immovable atoms positioned as a continuation of 
the crystal. This procedure also gives the bombarded face some character- 
istics similar to a binding energy. 

In order to simplify the calculations, we defined a quantity called 
the "lattice unit". One lattice unit is equal to half the length of a 
cube side. For copper, the length of a cube side is 3.614A, therefore 
one lattice unit is 1.807A. A three dimensional Cartesian coordinate 
system is used for reference with the origin placed with a (100) face of 
the movable 63 atom “core” parallel to the x-z plane at y=0. Any penetra- 
tion into the lattice core from this plane is in the + y direction. The 
immovable atoms surrounding the core are in the x-z plane at y=-1 and y=5, 
the x-y plane at z=0 and z=6, and the y-z plane at x=0 and x=6. For conven- 
ience all atoms are numbered. Number one is the “bullet”, numbers 2-64 are 
the "core" and 65-172 are the immovable atoms. (see Fig. 1) 

The program depends upon a potential function subroutine; so potential 
modification is accomplished without any changes in the main program, The 
Born-Mayer and the Thomas-Fermi-Firsov have been used in the present work 
(see Appendix VI for exact forms and constants). All calculations are 
with the potential in "eroded" form. To "erode" a potential, the value of 
the potential at the nearest neighbor separation is subtracted from the 


value of the potential at smaller distances, The eroded form makes the 
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potential zero at distances greater than nearest neighbors. 

The interaction between atoms is represented by an eroded repulsive 
potential. Also, the entire lattice is initially at absolute zero since 
no vibrational energy is simulated. These two approximations imply that 
the lattice has no potential or kinetic energy before interaction with the 
bullet and all the energy in the lattice at any time thereafter is derived 
from the bombarding atom. 

We realize that each atom in the lattice is contained in some type of 
“potential well", but we are measuring the energy from the ground state 
level of this "well" rather than from the true zero of potential at infinite 
separations, Thus, before interactions with the bullet, our lattice is not 
held together by the surrounding stationary atoms as in the Gibson et. alee 
model because the erosion also has removed all forces. Forces appear only 
when atoms move from their equilibrium positions. 

The use of an eroded potential does not provide as good a model as the 
one used by elisa et. al., but it reduces the computer time required for 
calculations. 

The "binding energy" created by the layer of fixed atoms covering the 


front face is probably much too large, but it does approximately simulate 


the behavior of the next layer of atoms. 


Newton's equation of motion can be rearranged to give the change in velocity 
of a body acted on by an unbalanced force. ( i AT Yn = L\/ ) 
The change in velocity can be related to a change in position if an average 
velocity is assumed (Xew =[V+ S¥]-AT + Maisie The unbalanced force 
used in the first equation is an average force calculated by a double itera- 


tion procedure as follows: (1) assume an atom at position 1 with velocity 1 


The calculations which move atoms are similar to those of Gibson et. 4.3 
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(2) calculate the total force on the atom as a result of all the other atoms 
in the lattice (this means normally only about 8-10 nearest atoms because 

the potential is eroded) (3) call this calculated force, force 1, and use the 
equation of motion to move the atom to a temporary position, position 2, 

(4) now repeat the force calculations for position 2, call this force 2, 

(5) go back to position 1, and use the average of force 1 and force 2 to 

move the atom to a new position, position 3. Procedures 1 through 5 consti- 
tute one "time step". Forces are eroded in the same manner as potentials 

and are calculated by a subroutine based on the partial derivative of the 
potential function with respect to distance, 

The basic value of AT is dependent upon the original energy of the 
bombarding atom. AT is the time in seconds required for the incoming atom 
to traverse one lattice unit. The atom will lose energy as it interacts 
with the lattice but MT remains unchanged. To add flexibility to the 
program, the basic value of AT can be multiplied by any desired factor 
(time step multiplier). AT is established at the beginning of the program 
and is constant for the duration of the calculations, 

The bombarding copper atom is originally located in the planey = al 2 
and given a velocity in the + y direction only. This requires the bombarding 
particle to approach perpendicular to the (100) face of the movable lattice 
core, The immovable atoms in the plane (y = -1) that cover the bombarded 
face do not interact with the bullet. 

An impact area was chosen en the (100) face so that all points in this 
area would be representative of any point in a (100) plane. An impact point 
is defined as the location on the (100) face toward which the bombarding atom 
is directed. It may or may not actually pass through this point. By moving 


in 1/10 lattice unit increments, 36 impact points in the impact "triangle" 








were assigned. (see Fig. 2 and 3). 

For each case the two body problem is solved after the many body 
problem with the same double iteration scheme, Lattice atom number eight 
was assigned as the "target"’ in the many body lattice problem so that its 
motion could be compared with the simple two body problem. The geometrical 
relationship and physical constants of the two body problem correspond to 
those of the incoming atom and "target" in the crystal. 

The program was written in FORTRAN language for use on a CDC 1604 
computer. The input consists of the incoming particle mass and energy and 
the time step multiplier. From these inputs the bullet is assigned a loca- 
tion and an appropriate velocity component in the + y direction. Also cal- 
culated is the actual value of AT. 

The next section of the program assigns coordinates to all atoms in 
the crystal, movable and immovable, gives each atom a number and sets their 
velocities equal to zero. If the distance between any two atoms is less than 
J2 lattice units, the force between these two atoms is calculated and stored 
in memory. This is done for each atom and the summation of all the forces 
calculated gives the resultant force on each particle. 

The change in the velocities of all movable atoms as a result of un- 
balanced forces is calculated; then the original velocities of all particles 
and the calculated change in velocity are used to move atoms to temporary 
locations. The calculations to obtain resultant forces are repeated for the 
temporary locations, old and new forces are averaged and used to move all 
atoms to new locations. Forces and velocities are calculated as vector 
components using the coordinate system established for the lattice, 

After each movement of all the atoms, the kinetic energy of the tar- 


get is calculated and compared to its former value. The procedure is repeated 
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until the kinetic energy of the target reaches a maximum. The potential 

and kinetic energy of each movable atom is then calculated. The total 
energy of the entire lattice including bullet is also found and compared to 
the original "bullet" energy. This comparison gives a reasonable measure of 
the accuracy achieved, The two body problem is solved in the next section 
using the same iteration technique. Appendices I-V contain a detailed ex- 
planation of the program. 

For each "run", the computer will solve the many body and two body 
problem 36 times, corresponding to the 36 different impact points. The form 
of the printed output for each impact point is shown in Fig. 4 and 5. The 
"triangular" output as seen in Fig. 6 indicates how two chosen parameters 
vary over the impact area, 

The computer running time for each of the 36 impact points was 2-4 
minutes, the time varying with the value chosen for AT. We originally 
Assumed that the accuracy of the results would increase for smaller and 
smaller values of AT. We found that the value of AT which will give the 
minimum error in total energy appears to be a complicated function of the time 
step multiplier, the incoming bullet energy, and the impact point. By trial 
and error methods, we found that below 100 ev for the Born-Mayer potential, 
a time step multiplier of around 0.09 would produce reasonable errors at all 


impact points. 


11 





4, Potentials 

The Born-Mayer potential is an exponential function of the internuclear 
separation with constants which may be obtained from experimentally measur- 
able elastic moduli. The constants used in the program are the same as 
those of Gibson et. aie with Edicial emphasis on their potential number 
two. Near equilibrium separations, the Born-Mayer potential is thought to 
be an adequate approximation, We expect the potential to fail for energies 
above a few hundred electron volts. 

The Thomas-Fermi-Firsov (TFF) potential is the result of theoretical 
work by Firsov’’ based on the Thomas-Fermi model of the atom. It is a 
screened Coulomb potential with a more complicated screening function than 
that used in the Bohr potential. The TFF is satisfactory for intermediate 
separations. This places it tetween the Bohr potential used for small separa- 
tions and the Born-Mayer potential. Abrahamson and Hatcher? state that the 
Thomas-Fermi approximation becomes unreliable when the internuclear distance 
exceeds AVIA. The form of the TFF potential and the appropriate constants 
used in the program are the same as those used in computer programs at the 
Oak Ridge National Laboratory:’. 


The mathematical forms of the potential and force functions as they ap- 


pear in the program are derived in Appendix VI. 


wae reference 15, p. 1233. 


ao B. Firsov, J. Exptl. Theoret. Phys. (U.S.S.R.) 32, 1464 (1957). 
translation: Soviet Phys. -JETP 5, 1192 (1957). 


To a3 A. Abrahamson and R. D. Hatcher, Phys. Rev. 121, 159 (1961). 


oe T. Robinson (private communication). The author would like to 
thank Dr. Robinson for the use of his results prior to their publication. 
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5. Results 
A. General description of events in the lattice. 

The number of atoms disturbed in the lattice is a function of two 
variables, the number of time steps the computer performs and the impact 
point of the bullet. In general, the lattice is more disrupted by impact 
points in the area near x = 3.0, z = 4.0, and the lattice is least af- 
fected by head on collisions with the target. No more than one half of the 
atoms were ever disturbed before the target reached a maximum in kinetic 
energy. 

The impact area is not symmetrically located with respect to the sur- 
rounding immovable atoms but for appropriate points in the impact area, 
excellent symmetry of displacement was observed for the movable atoms. This 
indicates that the size of the model is adequate for present purposes. 

Because of the limitation placed on the number of time steps, only a 
few atoms have any significant kinetic energy when the interaction is stop- 
ped. These are normally the target, bullet, and the atoms directly behind 
the target and bullet in the lattice. 

Although our primary purpose was to compare lattice interactions with 
the two body interaction, we altered the program so that it would run a 
predetermined number of time steps in order to observe “chains” and “chan- 
nels", Impact points at x = 3.0, z = 4.0 and x = 3.0, z = 3.0 correspond 
to the beginning of (100)chains. These are not "close packed" chains but 
they transfer energy reasonably well, especially above 100 ev. Impact point 
x = 2.5, z = 3.5 is the beginning of a (100) channel. The (100) channel is 
not as "wide aeenias a (110) channel and a 100 ev bullet does not travel 


more than “#4 3A into the lattice. 
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B. Comparison of lattice interactions with the two body interaction. 

A strict comparison of interactions in the lattice to the two body 
problem is not possible, but certain limited correlations can be made. The 
difficulty arises because the end of the interaction in the lattice cannot 
be defined explicitly. There are several criteria available that may be 
used in an attempt to compare the interactions. If the two body interaction 
was allowed to continue until the potential between atoms was zero, for 
comparison purposes the lattice interaction could procede until: 1) the 
potential between target and bullet is equal to zero 2) the total energy 
(kinetic plus potential) of the target reaches a maximum or 3) the kinetic 
energy of the target reaches a maximum. All of the criteria apply tc the 
same physical situation in the two body problem but are not equivalent in 
the lattice interaction. 

We chose the maximum kinetic energy of the target as the comparison 
stopping point because this represents the time when the target and bullet 
have almost ceased to interact and are just beginning to interact with other 
lattice atoms. 

A fourth criteria was suggested after the present results had been 
obtained. When the kinetic energy of the target in the lattice reaches a 
maximum, find the distance between the target and bullet and allow the two 
body interaction to proceed until the distance between atoms peregunil to the 
separation in the lattice. This criteria is possibly better than the one 
chosen but at the present time it has not been utilized. 

The geometrical relationships of the target, bullet, and lattice intro- 
duce complications into the comparison scheme. The motion of the target in 
the lattice is not restricted isotropically (see Fig. 7). If the inter- 


action with the bullet moves the target toward atom number six (shaded area), 
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the movement of the target is restricted more than if the target motion is 
initially toward the cross-hatched area. Restriction of the target motion 
lengthens the interaction time between bullet and target, therefore the 
energy transferred to the target is greater than in the two body problem. 
This effect is noticable only at impact parameters of 0.5 lattice units or 
greater, At impact distances smaller than this the interaction occurs so 
rapidly that the target does not move a significant amount and the effect is 
reduced, The kinetic energy transferred at distances smaller than 0.5 lattice 
units is less than that transferred in the two body case because there is 
absorption of energy by the remainder of the lattice. The geometry of the 
bombarded face is not the only factor that effects energy transfer. The 
atoms. behind the target in the lattice also play an important role but their 
effect is not as immediately obvious as those just discussed. 

Figure 8 illustrates some of the points just mentioned. Along the line 
AB (x = 3.0, z = 3.0-4.0) the target behavior in the lattice is similar to 
the two body interaction. These impact points cause the target to move in 
the negative z direction toward a "hole" in the lattice face. At point D on 
the line AC , the effect of restricted target motion becomes obvious, and 
also in the area enclosed by the dashed line. The geometrical effects of 
the lattice are fairly predictable and no unusual phenomena are observed. 

For a static system of two equal mass atoms with a conservative repul- 
sive force between them, the total potential energy of the system can be 
halved and the result assigned as the energy of the atom. If the two atoms 
are allowed to move, the total energy (potential plus kinetic) of each will 
be constant and equal to the original potential energy. In a system composed 
of three or more atoms, this is generally not true. For any system where the 


atoms are not static, the process of assigning an energy to any one atom is 
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no longer possible (there is one exception). This is the situation in the 
lattice. 

The potential energy of an atom in the lattice is defined in the program 
as half the potential associated with its position. This definition is merely 
a convenience. 

The area of the impact triangle from point A to the arc E-F (see Fig. 8) 
is least effected by the geometry of the lattice and closer inspection of 
these interactions is justified. The kinetic energy transferred to the 
n-body target is less than that of the two body target. In percentages, 
the kinetic energy transferred becomes less as the bullet energy is decreased. 
(see Figs 8 through 11). Im each case the target also acquires some poten- 
tial energy (see Figs. 12 through 15). As stated earlier, it is difficult to 
assign a specific potential energy to the target. We do know the limits 
involved, i.e. the target could eventually receive none of it, or twice the 
indicated value. Either limiting situation is unlikely. 

If we make the assumption that one-half of the potential "belongs" 
to the target, then the total energy of the target is approximately equal to 
the energy transferred in the two body problem (within 2%) for energies great- 
er than 50 ev. 

Bullet behavior for the area under discussion is also very similar to 
the two body problem (see Figs. 16 through 19). For the remainder of the 
impact triangle, the bullet is strongly effected by other atoms in the lat- 
tice, 

Three runs were made with the TFF potential at energies of 5, 10, and 
30 kev (Figs. 20 through 25). Since the TFF potential is unreliable beyond 
1A, the section ABC in Figs. 20 through 25 is the only part of the impact 


area that can be considered. This section includes all impact points that 
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result in a closest point of approach of 1A or less, We noted the same 
general agreement with the two body interaction as found for the Born-Mayer 
potential. The only apparent difference is a shift in the energy scale, 

For a head-on collision in the lattice there is very little geometrical 
effect, and an "effective" mass can be described. At 25 ev, the mass of the 
target is apparently about 2.5 times the bullet mass, but the energy trans- 
fer is still approximately 80%. For a 500 ev bullet, the effective mass 
increase is negligible. We used the program, as modified to observe chains 
and channels, and could find no evidence of bullet recoil which could be 
attributed to an effectively heavy target as proposed by Henschke’, Although 
an effective mass concept for the head-on collision is possible, it does not 
properly describe the subsequent motion of either target or bullet. If the 
bullet motion results in a glancing hit with the target, no single effective 
mass can be assigned to the target because it is a function of the impact 
parameter. An average effective mass is not applicable because the geometri- 
cal effects of the lattice on the direction and energy of the recoil atom are 
far more significant than the mass of the target in a glancing hit. 

These results apply to a collision anywhere in the lattice and are not 
limited to the interaction of surface atoms with incoming particles, The 
results indicate that fenechke’ se: apparent mass concept and the rebound 
phenomena associated with it are not a good description of collision events 


in a lattice. 
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6. Conclusions 
A. The Bullet. 

For energies above 40 ev in the impact area AEF (see Figs 16 and 17), 
the n-body bullet's kinetic energy after interaction is essentually the 
same as the two body bullet. Between 40 and 400 ev for impact points out- 
side the area AEF, the n-body bullet has considerably less energy after inter- 
action than the two body bullet. Above 400 ev, the n-body and two body 
energy transfers agree within 3% for all points in the impact area. 

The angular behavior of a bullet in the lattice is much more complicated. 
Even at energies above 400 ev, the scattering angle in the n-body problem is 
affected by the geometry of the lattice. In general, the scattering angle in 
the lattice is smaller than that of the two body problem. [It is possible that 
the lattice is attempting to focus the energy into preferred directions. Our 
bullets in the (100) direction appear to focus in the (100) direction, especial- 
ly at higher energies. 

The cummulative effects of energy transfer and lattice geometry will 
affect the range of energetic atoms. The reduced scattering angle in the 
n-body model for atoms with energies above 400 ev should lead to ranges that 
exceed those found by two body approximation methods. Below 400 ev, the 
bullet loses more energy than in the two body case; so the effects may cancel, 


or perhaps ranges calculated by two body methods may be too large. 


B. The Target. 

In the area AEF (see Figs. 8, 9, 12,13), the n-body energy transferred 
to the target is the same (within 27%) as that in the two body problem for 
energies above 50 ev. For the remainder of the impact area, the n-body 


target receives more energy than the two body target for all energies above 
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25 ev. Isolated points at various energies may not conform to this rule but 
the effect is always present for impact points near the apex of the impact 
area (x=2.5, z=3.5). Near the apex, if the bullet energy is large (300 ev 

or greater), the energy transferred to the target is a very small fraction 

of the total for both the n-body and two body problems, However, the energy 
transferred to the target in the two cases can differ by as much as 75% of the 
transferred energy. 

As a consequence of this phenomena, the low energy portion of the n- 
body target atom energy distribution (assuming more than one collision has 
occurred) will be much higher than that expected with a two body collision 
assumption for lattice interactions. Precise measurements of energy trans- 
fer for large impact parameters do not appear to be warranted in view of the 
results obtained. 

Unfortunately, the face centered cubic lattice structure was incorporat- 
ed in the program before we learned of vekater’ sce experimental work with 
molybdenum targets (bcc), but certain qualitative comparisons are still 
possible. Veksler has interpreted the lattice behavior in terms of the ef- 
fective mass concept but his general conclusions are consistent with our 
work, Our model gives good evidence to support Veksler's’ position that the 
pair collision model using elastic spheres is not acceptable. 

Two body approximations inherently imply that certain information is 
available from preliminary n-body calculations. This work attempts to answer 
certain questions about interactions and cross sections. As anticipated, the 
results are not definitive but they do indicate sensitive areas which require 


further examination. 


oe I. Veksler, J. Exptl. Theoret. Phys. (U.S.S.R.) 42, 325 (1962). 
translation: Soviet Physics - JETP 15, 222 (1962). 
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Fig. 1. THE MODEL 
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Fig. 2, FACE (100) of LATTICE CORE. 
Numbers indicate atoms; shaded portion is impact area 
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Fig. 3. IMPACT AREA 
Numbers indicate impact points. 
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Fig. 4. TYPICAL COMPUTER OUTPUT for an impact point 
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Target movement is more restricted if initial motion is in shaded area 
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APPENDIX I 
GENERAL DISCUSSION OF THE PROGRAM 

The masses of the lattice atoms are equal to that of copper and are 
an integral part of the program. The bullet mass however, can be varied by 
changing the input data. The potentials and forces are calculated by func- 
tion subroutines at the end of the program. The main program is not dis- 
rupted when a different potential function is substituted. 

A coordinate system was established with an atom located at the origin. 
For a face centered cubic lattice in this coordinate system, the sum of the 
coordinates of an atom is always an even number, Fixed point numbers were 
used for lattice units, and the volume which contains all the atoms (movable 
and fixed) was scanned in the x, y, and z directions. If the sum of the 
coordinates of a point was even, an atom was placed there, otherwise a space 
was left. Truncation of quotients in fixed point arithmetic allows a simple 
test for odd or even numbers. The fixed point coordinates of the atoms were 
stored in memory as floating point numbers for later calculations. The 
"core" atoms were assigned numbers from 2-64 and the fixed atoms were number- 
ed from 65 to 172, 

In the calculation of the forces or potentials associated with an atom, 
the atoms closer than nearest neighbors in an undisturbed lattice must be 
considered. A "nested' Do loop was used for these calculations. The out- 
side Do loop was indexed from 1 to 64 and the inner Do loops initial index 
was always one greater than the outer Do and ended at 172. This arrange- 
ment avoids the calculation of forces or potentials twice for the same pair 
of atoms. The x, y and z coordinate separations, in that order, between 
the atom under consideration (index of outer loop) and any other atom were 


tested. If any of these separations were found to be greater than ROE, 
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the inner Do index was advanced by one. If the coordinates passed these 
three tests, the square of the distance between atoms was compared to the 
square of ROE (ROE2). The true distance of any atom that passed this test 
was used in the force or potential function and the result applied to both 
atoms. (In potential calculations, half the potential was given to each 
atom). If the square of the distance between atoms was equal to or greater 
than ROE2, the inner Do index was advanced by one. Most of the atoms will 
be eliminated in the first three tests. The procedure avoids time consum- 
ing square root calculations by the computer. 

The double iteration method used to move atoms requires two calcula- 
tions of resultant forces for each "time step". The same section of the program 
was used for both force calculations, but the equations required to move the 
atoms twice in the time step are not identical. The fixed point variable 
INDEX was used as a switching device to jump over the motion equations that 
locate the atoms at temporary positions. The jump is executed after the 
second calculation of forces, 

The movement of atoms is continued until the kinetic energy of the 
target is a maximum. [It is necessary to advance one time step beyond the 
maximum to establish the maximum. The particle locations and velocities at 
the time of the maximum are still in memory, and these are used in kinetic 
and potential energy calculations or as output. 

Although the two body problem is in two dimensions, the impact para- 
meter is the same as that in the lattice. The two body problem is solved 
in the same manner as the lattice except that the interaction is stopped 
when the atoms are separated by a distance equal to ROE which means the 
potential and force between atoms is zero. The interaction is initially 


Started with the separation of atoms slightly less than ROE. The AT for 
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the two body problem is much smaller than the one used for lattice calcula- 
tions. The small AT consistantly produces very much more accurate re- 
sults in the two body problem than in the lattice. 

The input consists of the bullet energy and mass, its x coordinate, 
and the time step multiplier. In order to calculate data for all 36 impact 
points, six input cards are required; each is identical except for the x 
coordinate of the bullet. The x coordinates assigned are 3.0, 2.9, 2.8, 
2.7, 2.6, and 2.5 in that order. After receiving the input from one card, 
the corresponding z coordinate on the line AB is calculated (see Fig. 6). 
The computer proceeds through the program for this point, then subtracts 0.1 
lattice units from the z coordinate and repeats. It procedes down a column 
until all calculations for impact points on the line BC have been made. The 
computer pauses at this point to await instructions. Here it may receive a 
new data card and start another column. After the calculations for the point 
on the apex of the triangle (x#2.5, z=23.5) are complete, the computer can 


be directed to print the output in the triangular form shown by Fig. 6. 
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APPENDIX II 


BLOCK DIAGRAM 


Dimension and common 3 
Statements 


wBroded force and potential 
subtraction terms =, 
BZ and stopping point 


Original location of bullet ie 


Radius of effect (ROR) 
—Cinput _) 





BLOCK A 


Velocity of bullet 
Time step 
Constants used in velocity 
and movement equations 
Sst lattics velocities =0 








BLOCK B 


Location of lattice 
BLOCK C 


set potential and kinetic energy 
of lattice atoms=0 
Store original voositions offal 


atoms 
BLOCK I Dp 


To BLOCK & 


48 













From BLOCK D 






Set all forces = 0 
Calculate resultant force 
on all atoms 






Movement of atoms 
to temporary 
positions 


BLOCK Fl 





Movement of atoms to 
new positions 








BLOCK F2 






Test 
inetic energy 
of target for 
maximum 
BLOCK G 








Calculation of individual 
and total kinetic and 


otential energy 
, BLOCK H 


Output | 
See Fig. 3 prock J 
To BLOCK. am 
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Prom BEuCEK J 


Location of atoms for two 

body problem 

Calculate velocity of 
bullet 

set force =0 









BLOCK K 


















Test 
distance betwes 
atoms 
BLOCK Ll 








Calculate forces 
between atoms 


BLOCK L2 


Move atoms to 
temporary 
locations 


BLOCK Ml 


Move atoms to 
new locations 


BLOCK M2 
Calculation of ratios and 
parameters to comvars 


Outout two body vroblem with 
a lattice problem 
See Fig. 4 


BLOCK N 


To BLOCK R 
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From BLOCK N 


Storage of comparison 
data for each impact 
point 





BLOCK R 









Test 
for recycle 
BLOCK ASE 










Output 
Impact area in 
triangular and tabular 


form BLOCK W 





Bnd 


Potential 
function 


BLOCK Y 













Force 
function 


BLOCK Z 
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APPENDIX Il 


DETAILED FLOW CHART 





RoE= V2 
ROEZ=+ 2.0 


EPST= FoRC(RcE) 3 one terms subtracted from 
ic EUNemoN calculated potential and 


BLOCK A 


save memory locations for subscripted 
variables 


impact point number 


lattice unit distance between nearest neighbors 
and the square of that distance 


force that "erode'' them 


aA \Al 


3 INP LT ie starts here for every "column" 
in impact area 
EV, Bx, OTI, PMPSS in imp 
peo BX+ 1.0 


B2S= ico Rx 









RKC = BX 
RZ (1) = Be 
PKEGR=-1000 
CYCLE = O.O 
INDEX= O 


RY(I)= -N2 


5 







TO 
BLOCK B 


starts here for every impact 
point; assigns coordinates 
to bullet 


2 





VKXCi)= O&O 


v¥(i) 2A [evaes54 45,51\x10) 
(000.0% PMASS 


V2) = 0.0 


pr = DTK) * 1.807 x10" 
RN= \05.4b3411 x \o7’ 
Dr coc = DT /3.b14 x io” 
DTORM = PT/RM 
DTO2RM = BTORM/2,0 
Biy| = PMASS * 1.6598 x 10° 
PTORMI = DOT/RM\ 
DTO2ZRM| = PTORMI ZC 


“> 









DO 250 1 2, 


VXE) = 0.0 
VY (1) = 0,0 
Vz(1)= 0.0 





5O 
CONTINUE 


TO 
BLOCK C 


a3 


BLOCK B 


velocity components of 
bullet 


calculates constants used 
in force and movement 
equations 


sets velocity of lattice atoms 
to zero 





BLOCK C 


M=2 first numbers to be assigned to lattice atoms 
Nats 2 for movable core, 65 for fixed atoms 
Te -| 
IT, JT, KT are X, Y, Z coordinates 
in the system used for reference 
DO bO J=1,7 k-(60) 


KT=O 


each point with integer values for 

coordinates is tested; if the sum of 
the coordinates is an even number, 
an atom is assigned to that position 


DOS4 Ke 1,7 


 - 


© 


IT=O 


B) 


DO58 I= 1,7 


NTC= IT+IT + KT 


QO 
Q 
oO 


% 
* 
> 

K 

0 


these decision processes 
determine if an atom is in the 
"core"; this is necessary in 
order to have core atoms with 
consecutive numbers 2-64 


A 
O 


V 
0 


Lonlokin 
© ® ®® @ 











R2(N)= KT 











RX(M) = 17 
RY (M)= IT 
RZ(mM) = KT 


TT=IT+\ 
CONT | NUE 
CT = KT + 
54 


CONTINUE 


TT2IT+ | 


lo O 
end of Do loop over CONTINUE 


Y coordinates 


30 
BC Nie Ll 
RY (N)= UT 


BLOCK C ( con't) 


for a given Y and Z, X 
coordinates are checked 
first; then Z is changed and 
X coordinates checked again; 
after tests on X and Z, 

Y coordinate is changed 


end of Do loop over X coordinates 


end of Do loop over Z coordinates 


TO 
BLOCK D 





DO 6S T=&S,112 
PPE (1) = 0.0 

Rix (1)= RX) 
RLY (1) RY (1) 
Riz(I) = R2 (1) 








lb 
CONTINUE 


po bbl T= |, bH 
PPE (I) = 0.0 
PKE (I) = 0.0 
RXR (I) = RX(1} 
RYR (I) = RY C) 







RZR(1)= Re(1) 


BLOCK D 


sets temporary and permanent 
coordinates of immovable atoms 
equal to eachother 


Sets initial potential and kinetic energy 
of core atoms equal to zero 


stores original positions of bullet and 
core atoms 





BLOCK E 





70 

De cised=',64 
FX(I) = 0,0 
FY (1) = 0.0 





initial forces between core atoms 
set to zero 








1S 
CONTINUE 


DO 200 T=1\,b4 


TPi=1+\ 


Do200 J=IPFI,N 
keeps bullet from interacting 
with immovable atoms covering 
S-> bombarded face 


zo 
= DRX= RX(J) -RX(T) 
= tests X coordinate 
(200) ADRK = AGSE (DRX) separation between atom 


I and J 
© 
ADRX- ROE G00 | 
|<O 





IV 
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LO BLOCK E (con't) 


PRY= RY(J)-RY(I) 


n 


>o test Y coordinate separation 
ADRY —- ROE >— between atom I and J 
| <O 


DR2 = Re) -RZ(1) 


ADRZ= ABSF (DRZ 


S test Z coordinate separation 
= © 
ADR2-ROE = Geo) between atom I and J 


£0 
DIST= DRx7+ DRY + DRE* 


test square of distance between 
20 atoms; if square passes test, 
actual distmce is found and used 
| to calculate force 
<O 
DIST = / \DIST 
| FORCE 
FORCE® FORC(BLST)-EPST ke Fev wet ion 


FOD = FORCE/DIST 
FA = DRX»« FOD 
FB=DRY » FOD 
cc =DRZ * FOD 





reduce total force into X, Y, and Z 
components 









BLOCK E (con't) 












FX(J)=FXG)+FA 
FEY(J)=FY(U)+FB 
F2(J)=F2(T)+FC 







add force to total forces on atoms 
I and J if they are in the core 


Ex(1) = FX(I) - FA 
FY (1) = FY(I) -FB 
Emi ezia) ~Fe 


— 200 
end of Do loop for finding forces 
CONTINUE 


Do 23S I= 164 






round off procedures in 
computer may lead to very 
small resultant forces in 
situations where the forces 
should be zero; this procedure 
sets to zero any resultant force 
less than 10 ““ newtons 


23S TO 
CONTINUE BLOCK F4 


ag 





BLOCK F1 


DO 200—601> 2b 
RAX(I) = RX(T) 
FAX(1) = FX(I) 
Vix (I) = VX) 
a\so Same procedure 
For u & @ Components 





store position, velocity, and force 
of core atoms 







VX (I) = FX(I ) * DTORM + VX (I) calculate temporary velocities 
also same calc, for us 2 compone 


RXCT)= (VK(T) 4+ VLK(T))ADTOC + RIX Gq) move core atoms to temporary 
also Some ov 4 $2 coordinectes locations 


3200 
COW TINUE 


RIXQ)= RxXC) 
FAx(4) = FX(C4) 
Vax(d) = VX C4) 


also same procedure 


Lor 









store position, velocity, and force 
of bullet 







$ + components 
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BLOCK FI (con't) 


Vx(4) = EX(4) * DTORM + VX (4) 


Also same for Y b2 components 





since mass of bullet can be 
varied, equations of motion 
for the bullet must be 
calculated separately 







RX(4) = (VAC) + VEX(4))¥DTOC + RAX CL) 
also same for UY $2 coordinates 





INDEX= | 






GoTo 70 
(in BLOCK E) 





recycles to calculate forces at 
temporary locations 
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BLOCK F2 
FROM INDEX check 
im 


BLock Ft 


forces at temporary locations have 
been calculated 







new velocities calculated 
using average forces 


DO 425 T= 2 4 
VX) =[FXC) + F1xX(1))4DTO2]RM +V4X(I) 


RK (T) = (VX) + VAX(I)] 4 DTOC + RAX(Z) 
algo same caleulations Lor ue Z Components 


425 
CONTINVE 
NX (4) = [ExXG) + FAXG)J*PTOZRMI + VIK(1) | ce carate bullet 
RX (4) = L\ X(4) + VAX (4) | ¥DTOC + R1xX(1) calculations 


alse Same caleuletions Lov. Us? ComponentS 


INDEX=0 
CYCLE = CYCLE +\.0 


TO 
LOCK G6 





new positions assigned 
using average velocities 









BLOCK G 
- =@ 2 z 
PKES= 32,9 * 10 1Vx(8) fe VY (8) +\2(8)] kinetic energy of target 


NY 
oO 


PKESR- PKES kinetic energy of 
target compared to 
former value 


ZO 


PAE GR = PKES 





GeauTO aS 


if kinet 
(IN BLOCK E) 1 Netic energy is 


less than former 

: value, movement 
if kinetic energy of target greater than of atoms is stopped 
former value, recycle to calculate forces 

and move atoms again 

CONTINUE 


TO 
Brock H 
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BLOCK H 






DIX = RAX(@)- RIX) 
DLY =R1Y(8)-R1Y(4) 
prz =Riz(8)-R12(1) 







coordinate separation between target and bullet 





E e 2 2 
DBI =DIX +DIY + biz distance between bullet and target 
is found and potential energy calculated 


| POTENTIAL 
PoT1B=0.0) |PoT1S =PoTF(pal) “EPSTPL 7) conection 


PKEA2 0.518 * 10? x(V4x) + VY (1) + V2 (4) Js PMASS 


kinetic energy of bullet 









PKE(4)= PKE 
C= O1e 
TPKE=0.0 






DO b0O I= 1,&4 


calculate potential energy of each atom 
1Pi=1+1 and the total potential energy 


DO bOoO JI 2 


E> Fo) 


O 





SASRXE R1X(T)—RIX(I) 


ADRX= ABSF(DRx) 


BLOCK H (con't) 


prevents interaction of buLet with 
immovable atoms covering bombarded 


<O 


IV 


O 


DRY = RAY (J)-RIY (1) 


ADRY — ROE 


<O 


DRZ=R1IZ(5)- R4zZ(T) 


<O 


26 


NY 
O 


face 


test X coordinate separation 
between atom I and J 


test Y coordinate separation 
between atom I and J 


test Z coordinate separation 
between atom I and J 





BLOCK H (con't) 


<O 
DIST = DRX*+ DRY “4+ DR2* square of distance between atom 


I and J 

=2(0 | | 

= test square of distance, if it passes 
test, find actual distance and use it 
to calculate potential between I and J 





POTENTIAL 
POT=POTF(DIST) -EPSTP [7] cuverion 







PPE(I)= PPE(I)+0.5*POT half of potential between atoms is 
PPE (J)= PPE (3) +0,.5% POT assigned to each atom (I and J ) 





POT = TPOT + POT potential energy summation 


00 end of Do loop for calculating potential energy 
CONTINVE 


COL3= PPE(1)-O0.5 xPOT1IS 
PE BA= PPE (8)-0.5*PoT4i8 














DO 650 T=2,67 , ' : 
PKE (I)=32,9b x10 9(V1x(I) + VIY(I) +v4z(T) J 


TPKE = TPKE + PKE(T) 


TPKE = TPKE + PKEI total kinetic energy and total energy 
WS i= soem 


2 = oe 
PERE: | TE= Ev x |O percent error in total energy 


ZO 
Brock J 


kinetic energy of core 
atoms not including the 
bullet 
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BLOCK J 


OUT PUT 
see Fig. 4. 
Teo 
BLocKk K 


BLOCK K 


Conv = O.O\/DTTI change constants for two body problem; 


DTOC = DTCC * CONV DTI always equal to 0.01 in two body 
DTORM = DTOKM #CONV problem 


DTO2Z2RM =DT0O2RM« CONV 
DTORM 1 =DTORMIxX¥ CONV 


PTOZRM1=DTOZ2RMI*¥CONV 


RX2= 0.0 target located at origin 
RY2=0.90 
DSQ= (Bz -30 * + (BY- 3,0) square of impact parameter in the lattice 


RY1 = N DSO assign bullet Y coordinate to give proper 


impact parameter 


j—. 


XT P= RY 
~ =\RoE2-DSG =} calculate X coordinate of target to insure 
1="\Roz2- DS lO separation of bullet and target is slightly less 


than nearest neighbors separation 
VX22= 0.0 FA2= 0,0 

VY2= 0.0 FY1= 0,0 
VY1= 0,0 FYe = 0,0 
EX1= 0.0 


E\«.Ob354 4 ae 
Xie 5.51x\0 INDEX=O 
x ‘36 ASS. Biock U4 


assign velocity to bullet 
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BLOCK Ll 


\0O5 
DX = RX1- RXKe 
bY= RY1- RY2 


calculate distance between atoms 


3 is test to see if less than ROE 


DIST= DX* + DY? 


V 


O — 
BLOCK L2 


DIST = \ Das Tt 
FORCE calculate force between 
FORCE = FORC (DIST)- Bie ae Fuwcr ion atoms 


HA 


Fx4 = (Dx/ PIsT)¥ FORCE total force is divided into X and Y 
Exe “— FX] components 
FY1 = (DY/Dist)x FoRcE 
FY2:-FY] 
wed oO JBLock m2 
SO 
aa 
BLOcK M| 
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BLOCK M1 
RX1T = RX1 


RK2T = Ree 
EXIT = Fxi | store old positions, forces, and 
VXIT = VX] velocities 
VX2T = VAY 
also same for u 
Com ponents 








We? EXD SET ORM & VT calculate temporary 
VX¥22 FXK!N* DTORM +VXKZ velocities and positions 


RAL = LVx1 + VKIT)¥ DTOc + RKAT 
RX2:2 [VK2 *VK2T] % DTOC + RYOT 


also Y Components 







CO To recycles to calculate forces 
loos at temporary positions 


\woex = | 
Brock 1.1 









BLOCK M2 
FROM INDEX TEST - BLOCK L2 


V1 = (CEX1 + FXIT )x DTOZRMA+ VXLT 
NX2 = (FX2 + FX2T) eDTO2RM + VX2T 
RX1 = (Vx z VxaT ) AUTOR Haar average forces and velocities 
RAY = Cre +NX2T) * DTOC +RX2T are used 








calculate final velocities 
and positions 


also Y Components 


26 sa 
\NDEX =O \0OS recycles to distance test 
BLocK L141 
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BLOCK N 


TEST Bloce L] 





P2KE1 = 0,518x\08(VX17+ VY1*)xPMASS kinetic energy of target and 
P2 KE2=32. Alex 15° (vx 2% + V2?) bullet 


TKE = PekKEd + Pekea total energy of system 


Olimar TT 
See Fig. > 
CALCULATION OF SEVERAL 


VALUES AND RATIOS 
USED TO COMPRRE 
LATTICE WITH TWO BODY 


OUT PUT TO 
seas ri6. S BLOCK RK 









BLOCK R 










STORE VARIOUS VALUES 
FouNne FoR rattice AND 

TWO BODY PROBLEM IN TWO 
Dimensional ARRAY (C(TL)] 


index I is value stored 
index L is impact point number 


impact point number increased by one 
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: B2=B2- O.1 


>O 


SENSE 
SwitceH i 


DowN 


ao 
BLocKk UY | 






OUT PUT 


OF DATA 
see FiG, b 


IN 
O 


UP 


"TrrAn6uLar PRESENTATION 





BLOCK T 


Z coordina te of impact point lowered by 
0 10lattice units 
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Z coordinate of impact point 
tested to see if itis the end 
of a column 


if sense switch one is up, new 
data card is read to start another 


column in "triangle" 


BLOCK W 





: BLOCK Y 
FUNCTION POTF(DIST 


COMMON 


eae, AFAC, ae Gy 





Aes 


, -3.7% 
Pore = (ETFP)(1.0+ [bist ¥aracy 5°34 ) : 
iain 


END 


BLOCK Z 


FUNCTION FoRC(pTST ) 
COMMON 
ETEP, AF AC ETFF 


Y= (orsT» AraC) 
YP=1,0+Y 


ee 34 






Z=(YP 


orc = (er EF) «(Z)*(1,.0 + 2.9999 *YKp | 
Dist « DIST 


RETURN 
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APPENDIX IV 


DEFINITION OF VARIABLES 


RX(TI) X coordinate of atom number I at any time (lattice units) 


R1X(T) x coordinate of atom number I at former position (same as 
position 1 in procedure section), stored while average 
forces are computed (lattice units) 


RY(I), RZ(I), RLY(I), R1Z(I)  y and z coordinates defined in the same 
manner as RX(I) and R1X(I) 


FX(T) x component of force on atom number I at any time (newtons) 
F1LX(I) x component of force at former position of atom number I 
(newtons) 


feel), FZ(I), FIG), F1Z(T) y and z components of force defined in the 
Same manner as FX(I) and FIX(TI) 


VX(T) x component of velocity at any time of atom number I (m/sec) 
V1X(T) x component of velocity of atom number I at former position 
(m/sec) 


VY(I), VZ(I), VIY(I), V1Z(I)  y and z velocity components defined in the 
Same manner as VX(I) and V1X(I) 


C(I,L) matrix for storing data I from impact point L 
PPE (T) potential energy of atom number I (ev) 
PKE (TI) kinetic energy of atom number I (ev) 


RXR(I), RYR(I), RZR(I) x, y, and z coordinates of all atoms in initial 
positions (lattice units) 


ETFP, AFAC, ETFF constants used in Thomas-Fermi-Firsov potential 


L | impact point number 

ROE nearest neighbors separation (lattice units) 

ROE2 ROE squared 

EPST value subtracted from calculated force to erode the force 
(newtons ) 

EPSTP value subtracted from calculated potential to erode the 


potential (ev) 
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FORC 


POTF 


EV 
BX 
DTI 
PMASS 
BZ 


BZS 


PKE8R 


CYCLE 
INDEX 


DT 


RM1 


DTOC 


DTORM 


DTO2RM 


DTORM1 


DTO2RM1 


function subroutine that calculates force; argument is 
the separation of atoms (newtons) 


function subroutine that calculates potential; argument is 
the separation of atoms (ev) 


original kinetic energy of bullet (ev) 

original x coordinate of bullet (lattice units) 
time step multiplier 

mass of the bullet (amu) 

original z coordinate of bullet (lattice units) 


z coordinate of bullet at end of column in impact triangle, 
used as stopping point (lattice units) 


kinetic energy of target at the end of previous time step 


(ev) 

number of time steps 

switching device 

length of basic time step (seconds) 

mass of copper atom (kg) 

mass of bullet atom (kg) 

a constant used in movement equations; when multiplied by 

the sum of two velocities, it computes the average, multiplies 
it by DT and gives the result in lattice units (sec x lattice 


units/m) 


a constant used in velocity equations; it combines the 
quotient DT/RM into one constant to avoid repeated divisions 


(sec/kg) 

same as DTORM except it averages the force 

same as DTORM but used in bullet calculations only, since 
the bullet mass may be different from that of the lattice 


atoms 


same as DTO2RM but used in bullet calculations only since 
bullet mass may be different from that of the lattice atoms 


atom numbers 


y coordinate of an x-z plane 
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KT 
IT 
NTC 


DRX, DRY, DRZ 


z coordinate of an x-y plane 
x coordinate of a y-z plane 
sum of x,y, and z coordinates of a point in the lattice 


difference between x,y, and z coordinates respectively of 
any two atoms (lattice units) 


ADRX, ADRY, ADRZ absolute value of DRX, DRY, and DRZ 


DIST 


FORCE 


FOD 


FA, FB, FC 


PKE8 


DIX, DIY, DIZ 


D81 


POT18 
PRE1 
TPOT 
Tek 
POT 


COL3 


PE8A 


TE 
PERE 


CONV 


distance between two atoms or the square of the distance 
(lattice units) 


eroded force between two atoms (newtons) 


FORCE divided by the distance between two atoms (newtons/ 
lattice unit) 


x,y, and Z components respectively of the force between two 
atoms (newtons) 


kinetic energy of target atom in the lattice (ev) 


x,y, and z coordinate separation between target and bullet 
(lattice units) 


distance between target and bullet in lattice or the square 
of the distance (lattice units) 


potential between target and bullet (ev) 

kinetic energy of bullet at end of time step (ev) 
total potential energy of lattice (ev) 

total kinetic energy of lattice (ev) 

potential between any two atoms (ev) 


potential of the bullet with respect to the lattice minus 
the target (ev) 


potential energy of the target with respect to the lattice 
minus the bullet (ev) 


total energy of the lattice (ev) 
percent error in total energy (percent) 
time step multiplier conversion factor for two body problem; 


makes two body time step multiplier equal to 0.01 regardless 
of the value used for the lattice problem 
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DSQ 


XIP 


RX2, RY2 


RX2T, RY2T 


square of the impact parameter in the two body problem 
(lattice units squared) 


impact parameter in the two body problem (lattice units) 


x and y coordinates of target at any time in two body 
problem (lattice units) 


x and y coordinates of target at former position, stored 
while average force is computed (lattice units) 


RX1,RY1,RX1T,RY1T x and y coordinates of the bullet in the two body 


VX2, V¥2 


VX2T, VY2T 


problem defined in the same manner as RX2,RY2,RX2T, and 
RY2T 


x and y components of target velocity in two body problem 


(m/sec) 


x and y components of target velocity in two body problem 
at former position (m/sec) 


VX1,VY1,VX1T,VY1T x and y components of velocity of the bullet in the two 


FX1, FY1 


FX1T, FYIT 


body problem defined in the same manner as VX2,VY2,VX2T, 
and VY2T 


x and y components of force on the bullet in the two body 
problem (newtons) 


x and y components of force on the bullet in the two body 
problem at its former position (newtons) 


FX2, FY2, FX2T, FY2T x and y components of force on the target in the two 


DX 


DY 


P2KEl 


P2KE2 


TKE 
T2RA 


B2SA 


body problem defined in the same manner as FX1, FY1,FXIT, 
and FYILT 


x coordinate separation between atoms in ehe two body problem 
(lattice units) 


y coordinate separation between atoms in the two body problem 
(lattice units) 


kinetic energy of the bullet in the two body problem after 
interaction is complete (ev) 


kinetic energy of the target in the two body problem after 
interaction is complete (ev) 


total kinetic energy in the two body problem (ev) 
recoil angle in the two body problem (degrees) 


scattering angle in the two body problem (degrees) 
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SRE8& 
SRE1 
ATLRA 
TLRA 
ABLSA 
BLSA 
RBMT 


EABL 


resultant velocity of target in lattice (m/sec) 
resultant velocity of bullet in lattice (m/sec) 
quotient of two velocities; used to find TLRA 
recoil angle in the lattice (degrees) 

quotient of two velocities; used to find BLSA 
scattering angle in the lattice (degrees) 
bullet to target mass ratio 


energy absorbed by the lattice; does not include potential 
or kinetic energy of target and bullet (ev) 


ET2,ETL, RET, RKEP, RSA2L, RRAZ2L, RBFI ratios used to compare lattice interaction 


with two body problem; the ratios are in terms of variables 
already defined. 
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APPENDIX VI 


POTENTIALS 
1. Born-Mayer 


A. Potential 


b(r) = inf a ro) re 


r, = nearest neighbor separation at zero pressure and absolute zero 


i for copper = 2.51A 


The three different sets of constants used in the potential form are 


identical to those used by Gibson et. ai 
Potential A(ev fo 
1 0.0392 16.97 
2 0.051 13.00 
3 0.1004 10.34 


Potential No. 2 was used extensively in the program. It can be reduced 


to a simplified form in the following manner: 


-13.0(r-2ss\)/ 551 
D(r)= O,O3| «Cc 


if distances are expressed in lattice units 


(13,0 -9.208 r) 
QM (r)= OxO05 | x © 


in the potential function subroutine this becomes 
POTF = 0.051 * EXPF(13.0-DIST * 9,208) 
Potentials one and three become; 
1) POTF = 0.0392 * EXPF(16.97-DIST * 8.9996) 


3) POTF = 0.1004 * EXPF(10.34-DIST * 7.3115) 
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B. Force 


OLowd - pe rs ae 
or 


substitution of constants for Potential No. 2 


0 LO(r)] _ — (13.0)(0,08 1) a molr-zeslfassi 


or > 55 \ 





if distances are expressed in lattice units and electron volts changed to 


joules; 


a 13.0- 9.206r) 
Oo LOtr)] = — 0,4|/b14 « 10 x - 


y 


in the force function subroutine this becomes; 


FORC = 0.41614E-09 * EXPF(13.0-9.208*DIST) 


forces for potentials one and three become; 


It 


1) FORC = 0,3127E-09*EXPF(16.97-8.9996*DIST) 


2) FORC = 0.6507E-09*EXPF(10. 34-7. 3115*DIST) 


2. Thomas-Fermi-Firsov 


A. Potential 


vin = $e (3) 4) 


ot 





“he 
xX, = — - oe ZiE.C/ 


i 
an? \ 
Aa = 3,734 Ore == ( ¥ 0) Kay 
/: 
(% «EF 
les 
Oy - WK first Bohr radius of hydrogen 
= \\ 
= | — = 219233 ey Oren= AbOb2«l0 m 
let x = t/a. and ATFR = l/a,., (in lattice units) 


then x = ATFR * DIST (DIST in lattice units) 


(x) = f . (4 Jim = i - (ATER: ait) 


CX 
let ATFR&, = AFAC 


00,2034 shia 


on) = [lr (arc: ss 


V(r) = | | o(x) = Ere « OG) 
c XK AD: ATER. DIST 


= * 
let ETFP Ep! 2 * ATFR 


then 


34 


POTF = ETEP 


pros 


oles a 
2 + (AFAC * DIST) 
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B. Force 


set x = r/a 


EAC) _ Ww) dx 


Ox or 
ONO) — Ee}. 6%) , 1 d(om)], 
Sr US \? ie K | Ore 
Ka 
set Y = (% 
then 
dd (Oa) | = —%,%3 O(x) _Y_ 
d Xx (+ Y¥ 


O(VGr)) . Ere * 


- Ere *Orr || + 23 ee) b (x) 
or aan +! 


if the following substitutions and appropriate changes in units are made 


FORC = ETEE*z [1o+ 2.9999 
PIST « DIST : YP 
where ETFF = E 


sk 
tr * @qp /? 


x = DIST * AFAC (DIST in lattice units) 
0, G8 034 
= (x) YP= 129 ey 
~ 3134 
Z= (YP) 


ETFF and AFAC are defined at the beginning of the program - Y, YP and Z 


are calculated in the subroutine each time it is used 


a3 
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